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We compute non-perturbatively the evolution of the twist-2 operators corresponding to the average momentum 
of non-singlet quark densities. The calculation is based on a finite-size technique, using the Schrodinger Functional, 
in quenched QCD. We find that a careful choice of the boundary conditions, is essential, for such operators, to 
render possible the computation. As a by-product we apply the non-perturbatively computed renormalization 
constants to available data of bare matrix elements between nucleon states. 



1. INTRODUCTION 

The accurate knowledge of hadron parton den- 
sities is an essential ingredient for the experimen- 
tal test of QCD at accelerator energies. Their 
normalization is usually obtained from a fit to a 
set of reference experiments and is used for pre- 
dicting the behaviour of hard hadron processes 
in different energy regimes. The calculation of 
the normalization needs non-perturbative meth- 
ods. In order to have a phenomenological impact 
this determination must have a precision compa- 
rable with the experiments, and must have all the 
systematic uncertainties under control. In this 
proceedings we will mainly summarize the results 
obtained in pP , to which we refer for any unspec- 
ified notations. 



2. MOMENTS OF PARTON DISTRIBU- 
TION FUNCTION 

The moments of parton distribution functions 
(PDF) are related to matrix elements of leading 
twist r (r =dim-spin) operators of given spin, 
between hadron states h{p) 



{h{p)\0^,...^,\h{p)) = Af(^-i)0.)p^,.-. 
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(.(^-i))(M)=Af(^-i)(^ = Q) 
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On the lattice the 0(4) symmetry is broken to 
the hypercubic group H{4), and the 2 irreducible 
representations of the non-singlet operators for 
iV = 2 are 

r - lA - ir3 
044(2;) = 4>{x) [74 Di -3 Z.'^fc \ y^(2^) (3) 

k=l 



71 D2 +72 Di 



yV'(a;) (4) 



Our setup will be QCD in a finite space-time 
volume of size T x with T — L. We choose the 
same boundary conditions of |4l5j . namely inho- 
mogeneous Dirichlet boundary conditions at time 
xq — and xq = T and periodic spatial boundary 
conditions up to a phase for the fermion fields. 



iP{x + Lk) e"''=?/'(a;) 



(5) 



The strategy used to compute the non- 
perturbative evolution of the operators in eq. 
H3I4|) resembles the strategy used by the ALPHA 
collaboration to compute the running quark mass 
|S] . The evolution from initially large L (low fj) to 
small L (high ^jl) is obtained applying the so called 
step scaling function (SSF) (cfr sect. 3). Once the 
perturbative regime is reached (and this must be 
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Figure 1. Relative errors for the Z factor com- 
puted with 400 measurements on a 16^ lattice at 
g|^(L)=3.48 



Figure 2. Continuum extrapolation of the SSF 
for selected values of gsp{L) 



checked) one continues the evolution in pertur- 
bation theory computing the (scale and scheme 
independent) RGI matrix element. The connec- 
tion with experiments is obtained then, making 
the adequate perturbative evolution of the RGI 
matrix element in the MS scheme. 

3. RENORMALIZATION 

The renormalization conditions for the local 
operators are given by 

ORifi) = Z^\afi)0{a), Onifi = L^^) = O^") 

The correlation functions to compute the Z factor 
are 

fo{x,IL, 0) = -^ J2 P{x) §,(y, z)) (6) 



x,y,z 



u,v,y ,z 



(§,(y,z)Sl(u,v)) 



(7) 



where §^ and are suitable quark sources to 
probe the operators O. With this definition the 
renormalization constants are obtained by 



Z(a/ L, fJ.) — c 



fo{xo/L,0) 



fS\^o/L,9)' 



Optimal choice of 9 and xq is mandatory to ob- 
tain a reliable signal of the correlation function 
(jHJ . In fig. n] a study of the relative error of the 
Z factor is performed. A similar analysis can be 



performed for the cut-off effects and for the con- 
vergence of perturbation theory computing the 
2-loop anomalous dimensions for these operators 
in the SF scheme jj. From these studies a good 
choice turns out to be di = 1.0, 62.3 = and 
Xq = f . 

To map out the L dependence recursively we 
use the SSF, rigorously defined on the lattice by 



~ lim 'Lzoiu^a/ L) 



Szo(u,a/L) 



Zo{u,2L/a) 
Zo{u,L/a) 



(8) 



fspiL) (9) 



The values of [3 corresponding to a fixed running 
coupling are available in (5.. We have computed 
the SSF at 9 values of the renormalized coupling 
{gsp{L) = 0.8873 to 3.48) corresponding to a 
range of energies that are roughly between 300 
MeV and 100 GeV. In order to have a better con- 
trol on the continuum limit we have performed 
the computation with Wilson and Clover action, 
even if in both cases one expects 0{a) lattice arte- 
facts since the local operators are not improved. 
In fig. 121 the continuum limit of the SSF for some 
values of ggp{L) is shown. It is clear that a re- 
liable (constrained) linear extrapolation with a 
small slope is possible. 

The formula that summarizes the whole strat- 
egy is given by 



RGI 



hmJ^M^x 

a^O Zo{a,^lQ) 

<^Zo (m/MO , ^ (m) )^SF (ff ^ (m) ) 
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Figure 3. Continuum limit of the non- 
perturbative renormahzed first moment of the 
PDF in a proton 



where we use the n = 9 SSF computed with ji 

fJ-n 



— = cr( — ■ 
Mo Mo 



to jump from the non-perturbative scale fiQ to the 
perturbative (ultraviolet) scale /i. At this point 
one can try to do the perturbative matching using 



/ dx 


[7(2;) 
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computed with 3-loop [3 and 2-loop 7 functions. 
If the perturbative matching has been successful 
the quantity 



f^/w(Mo) = (j{^./^.a,g {^)) TsF{g (m)) 



(10) 



should be independent from the ultraviolet scale 
(ji. Indeed what we find is that the last 4-5 steps 
give a very nice plateaux (see fig. 10 in ref. P). 

So it is possible to continue the evolution from 
the last step using perturbation theory. 

4. PRELIMINARY RESULTS 

Using perturbation theory is possible to com- 
pute the Z^^^ , that is the fundamental quantity 
to relate bare matrix elements to any desirable 
scheme (e.g. MS) 
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5 Z.^G^(a) 



(0)j^(/i)%s(M) (11) 



(a) = Zo{a,^j.a) 
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It is then clear that knowing Z^^^ for a certain 
discretization allows to compute in our case the 
parton average momentum {x) in the proton by 



(x)(a) 



■oZ«G/(a)^_(^) 



(12) 



We then apply the Z^'^^ we have computed 
to the unpublished data |2] of QCDSF for the 
nucleon bare matrix element and to the published 
data (3j of LHPC available at only one value of 
(3. The continuum limit is shown in fig. |31 

5. CONCLUSIONS 

We have computed in the continuum and in a 
fully non-perturbative way the evolution of the 
twist 2 non-singlet operators with a very good 
precision (4%). This computation, combined 
with a calculation of the bare matrix element be- 
tween hadron states (for an application to pion 
matrix elements cfr. [H]), gives the renormal- 
ization group invariant matrix element (x) 
Then the RGI matrix element can be simply con- 
verted to any desirable scheme. The precision of 
the experimental data requires a better control on 
all the systematic uncertainties (non-perturbative 
renormalization, continuum limit, chiral extrap- 
olation, finite volume effects quenching). In 
this contribution we have shown how to have com- 
plete control over the non-perturbative renormal- 
ization and on the continuum limit. There is still 
a disagreement between the experiments and the 
lattice computation, but there are also still sys- 
tematic uncertainties in the lattice computation 
that must be carefully analyzed. It is clear that a 
comparison between experiment and theory can- 
not be reliably done with a lattice simulation at 
one value of the lattice spacing and without doing 
a non-perturbative renormalization. On the other 
systematic uncertainties works are in progress. 
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